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ABSTRACT 

With high-angular-resolution, near-infrared observations of the young stellar 
object T Tauri at the end of 2002, we show that, contrary to previous reports, 
none of the three infrared components of T Tau coincide with the compact ra- 
dio source that has apparently been ejected recently from the system (Loinard, 
Rodriguez, & Rodriguez 2003). The compact radio source and one of the three 
infrared objects, T Tau Sb, have distinct paths that depart from orbital or uni- 
form motion between 1997 and 2000, perhaps indicating that their interaction 
led to the ejection of the radio source. The path that T Tau Sb took between 
1997 and 2003 may indicate that this star is still bound to the presumably more 
massive southern component, T Tau Sa. The radio source is absent from our 
near-infrared images and must therefore be fainter than K = 10.2 (if located 
within 100 mas of T Tau Sb, as the radio data would imply), still consistent with 
an identity as a low-mass star or substellar object. 

Subject headings: astrometry — binaries: close — infrared: stars — instrumenta- 
tion: adaptive optics — stars: individual (T Tauri) — stars: pre-main sequence 

1. Introduction 

Stars frequently form in small, gravitationally-bound groups of three or more (Ghez, 
Neugebauer, & Matthews 1993; Mathieu 1994), usually situated so that each perturbs the 
motion of the others significantly. In such systems the stellar orbits develop chaotically. 
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Over time, gravitational interactions among the members and with the residue of the disks 
from which they formed, result in large orbital changes and even ejection, especially for the 
lighter components (see, e.g., Marchal 1990). As main-sequence stars are always seen in 
configurations that are quite different from those of young multiple stars, direct observation 
of orbital evolution in young multiple-star systems has long been sought. 

Thus the report by Loinard, Rodriguez, & Rodriguez (2003) was greeted with wide 
interest: one of the components of T Tauri has, in the last few years, undergone a dramatic 
change in orbital velocity, possibly resulting in ejection from the system. T Tau has long 
been the archetype of young solar-analog stars, and more recently has become an archetype 
of young multiple star systems. Since the discovery two decades ago that T Tau is not 
single, but has a companion 0'/7 to the south (Dyck, Simon, & Zuckerman 1982), high- 
resolution infrared observations using speckle interferometry or adaptive optics have been 
used to resolve the system into three stars and to measure accurately the motions arising 
from their orl)its: the visible, classic T Tau N, and two heavily-extinguished companions, T 
Tau Sa and T Tau Sb, separated by about O'.'l (Koresko 2000; Kohler, Kasper, & Herbst 
2000; Duchene, Ghez, & McCabe 2002). Meanwhile, radio astronomers used the VLA at 
2 cm wavelength to resolve the system also into two compact components, one coincident 
with T Tau N and a brighter one close to T Tau S (Schwartz, Simon, & Campbell 1986). 
The southern radio component exhibited motion that until recently was consistent with 
orbital motion around T Tau Sa (which does not appear in the radio images) and has been 
identified with T Tau Sb (Johnston et al. 2003; Loinard, Rodriguez, & Rodriguez 2003; Smith 
et al. 2003). The radio source is unresolved in VLBI images (diameter < 0.5 mas; Smith et 
al. 2003), variable in brightness, and exhibits strong and variable polarization in its radio 
emission (Smith et al. 2003; Johnston et al. 2003). This is the object that has apparently 
been ejected. 

Here we present new near-infrared observations which we use to show that the southern 
radio component cannot be the same as the star T Tau Sb. We suggest that T Tau has been 
a quadruple stellar system, with the lowest-mass member probably ejected during the past 
few years. 

2. Observations 

We observed T Tau under good observing conditions on 2002 December 24 with the 
Palomar Observatory Hale 200-inch telescope, the Palomar Adaptive Optics (PALAO) high- 
order image-correction system (Troy et al. 2000), and the Palomar High- Angular-Resolution 
Observer (PHARO) near-infrared camera/spectrometer (Hay ward et al. 2001). The brightest 
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visible component of the system, T Tau N (V = 9.6), served as the phase reference for the 
AO system. We also observed SAO 76481 as flux calibrator and point-spread-function (PSF) 
standard immediately after taking our images of the T Tau system. We chose the image 
scale to be 25 mas per pixel. 

In order to calibrate the plate scale and orientation of the array, we used our observa- 
tions of the multiple systems RW Aur (two components) and UX Tau (four components) 
and compared the measured separations and position angles with the corresponding mea- 
surements in White & Ghez (2001). We chose multiple systems whose components have 
separations wide enough that orbital motion should be negligible in the time period between 
the observations by White & Ghez (2001) (1996-1997) and ours (2002.98). The plate scale 
of our images comes out to 24.7 ± 0.4 mas per pixel, and the orientation of the detector 
array is such that north is up and east is to the left with ±1° accuracy. Thus position angles 
inferred from our data have uncertainty ±1°. 

Figure 1 is the image of T Tau which we took in a narrow-band filter (AA = 0.06 jim) 
in the infrared K band centered at 2.26 /im. (The components of T Tau are bright enough 
at near-infrared wavelengths that the use of the usual broadband filters is inconvenient, as 
will be seen below). A total of 20 frames, each with exposure time 1.8 seconds, were added 
to create this image. The resulting PSF is close to diffraction-limited, 102 mas in diameter 
(FWHM). The limiting 2.26 fim magnitude for parts of the image well separated from the 
brighter stars is 14.7 (5 a per pixel). Within 200 mas (two PSF widths) of the stars T 
Tau Sa and Sb, diffraction rings reduce the limiting magnitude to a range from 10.2 to 
11.2. All three of the previously-identified stellar components of T Tau appear in Figure 
1, well resolved from one another. Table 1 is a list of their positions and magnitudes. T 
Tau Sb lies 0'.'107, at position angle 289°, from T Tau Sa, and the magnitudes of the two 
southern components are about the same, with T Tau Sb brighter by a factor of 1.05. Since 
in late 2000 the corresponding separation and position angles amounted to 0'/092 and 267°, 
respectively (Duchene, Ghez, & McCabe 2002), we probably see orbital motion of T Tau Sb 
around Sa (see also Figure 3). 

The positions and fiuxes of the three observed infrared components of the T Tau system 
were determined by applying a PSF fit to the data. Standard IDL procedures were used to 
perform a Gaussian fit to the given PSF, the one of SAO 76481, and to apply it to the target 
stars. The position uncertainty resulting from this method is estimated to be 0.2 pixels, 
which corresponds to 5 mas. The fiux calibration was carried out by using the 2MASS K- 
band flux density of SAO 76481; we estimate the uncertainty of the absolute photometry to 
be ± 20 %, and that of the relative photometry to be a few percent. 

Figure 2 is our image (A = 2.145 yum, A A = 0.31 iim) of T Tau Sa and Sb, taken 



-4- 



5 minutes before our narrow-band K data. It is also a sum of 20 frames of 1.8 seconds 
exposure time each. Since T Tau N and our calibrator, SAO 76481, saturated even with this 
short integration time, we used the image of SAO 76481 observed with the 2.26 jjm filter to 
perform PSF-fitting on this image. The position difference between Sa and Sb is nearly the 
same in this image and the one taken in the 2.26 fim filter; the offsets listed in Table 1 are 
averages of the two. No absolute flux calibration was possible, but relative photometry was: 
T Tau Sb is brighter than T Tau Sa by the factor 1.36, somewhat higher than we obtained 
with the narrow-band data. In another image taken about 30 minutes later with the filter 
and an occulting mask over T Tau N, the fiux ratio of T Tau Sa and Sb was the same as 
in the previous image. The difference between this result and the flux ratio seen in the 
narrow-band K filter (Sb/Sa = 1.05) could be the result of differences in spectral features 
and degree of extinction between the two components. 

Comparing our narrow-band K data to the results from late 2000 by Duchene, Ghez, & 

McCabc (2002), T Tau Sa is fainter and T Tau Sb is brighter, each by about a magnitude 
at K. That T Tau Sb can be as bright as T Tau Sa in the K band means that the orbital 
motion of T Tau Sa relative to T Tau N is confused in all pre- 1997 measurements because 
of the presence of T Tau Sb. 

3. Analysis and Conclusions 

Figure 3 is a plot of the radio/infrared positions of the southern components of T Tau 
between 1983 and 2003, in a frame of reference in which T Tau Sa is at rest. The radio 
positions in this plot differ very slightly from those determined by Loinard, Rodriguez, & 
Rodriguez (2003). They are based on the radio/infrared registration of T Tau N and include 
the (very small) orbital motion of T Tau Sa with respect to T Tau N; the uncertainty in the 
distance between origins of the radio and infrared reference frames is about 15 mas (radial, 
RMS). The relative positions of T Tau Sa and Sb were measured directly in the infrared 
images. Our new position for T Tau Sb is quite inconsistent with anything along the track 
of the radio source between 1998 and 2001: it misses by 50-100 mas, much larger than the 
combined uncertainty of about 16 mas. For example, the linear extrapolation of the two 
most recent radio-source positions, the point indicated by a empty diamond in Figure 3, 
lies 78 mas from the 2002.98 observed position of the infrared star T Tau Sb. If a single 
object were to have followed the path described by both the radio and infrared positions, 
its velocity in the plane of the sky would have had to undergo major changes twice, once 
between 1995 and 1998, and once after 2001. Moreover, in at least one of the cases, the single 
object would have had no obvious partner with which to interact. Thus the probability is 
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much greater that the radio and infrared observations of the southern components of T Tau 
detect different objects. In anticipation of its detection in future infrared observations, we 
will refer to the radio source as T Tau Sc henceforth. 

Also appearing in Figure 3 is a model orbit, resulting from a minimum-x^ fit to the 
projected trajectory and transverse velocities of T Tau Sc measured at the VLA between 
1983 and 1996 (Loinard, Rodriguez, & Rodriguez 2003; Johnston ct al. 2003). In this exercise, 
we assumed T Tau Sa to be at one focus of the orbit. The fit to both trajectory and velocity 
improved after shifting all the VLA positions of T Tau Sc with respect to T Tau Sa by 4.5 mas 
to the east, an amount small compared to the uncertainty of the radio-infrared coordinate- 
system correlation. This small shift is present in the positions plotted in Figure 3. The 
agreement between model and observations is excellent; we suggest that the 4.5 mas shift 
is an improvement upon the previous determination of the infrared-radio reference-frame 
alignment. The parameters of the model orbit are as follows: period 19.1 years, semimajor 
axis 9.1 AU, assuming a distance of 140 pc (leading to a total mass of 2.IM0 for T Tau Sa 
and T Tau Sc). axis inclination 59 degrees from the line of sight, eccentricity 0.57, argument 
of periastron 168.5 degrees, position angle 0.5 degrees for the line of nodes, and periastron 
epoch 1997.39. T Tau Sb is a young Ml ± 1 star (Duchene, Ghez, & McCabe 2002), for 
which the mass is much smaller than 2.IM0; thus we have neglected it in this model. 

What is the nature of T Tau Sc? Its small size and large (kilogauss) magnetic field, 
indicated by the VLBI observations (Smith et al. 2003), are consistent with a star-like object. 
We did not detect any other object besides T Tau N, Sa, and Sb in our near-infrared images. 
At the position extrapolated for T Tau Sc at epoch 2002.98 (see Figure 3), our upper limit 
for the detection of a point source is about 10.2 magnitudes at 2.26 yum. Thus it is either 
fighter than the eighth-magnitude T Tau Sb, or more heavily-extinguished, or both. The 
velocities of T Tau Sb and Sc within the T Tau system are larger than those of T Tau Sa, 
which, if they have been gravitationally bound, indicates that both objects are less massive 
than T Tau Sa. Therefore T Tau Sc is probably an ordinary young star or brown dwarf, 
similar in mass or lighter than T Tau Sb, though quite different from this star in its magnetic 
field strength. 

The presence of two distinct objects, T Tau Sb and the radio source T Tau Sc, could 
provide a natural explanation for changes in each other's state of motion. As Loinard, 
Rodriguez, & Rodriguez (2003) point out, the path of T Tau Sc obeys Kepler's laws for an 
orbit around T Tau Sa from discovery through 1995, but departs strongly thereafter. The 
path of the infrared star T Tau Sb with respect to Sa shows a wide variation in area per 
unit time; between February 2000 and November 2000 the rate of area swept out is more 
than a factor of two greater than from late 1997 to early 2000, or from November 2000 to 
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December 2002. The curvature and orientation of this path probably indicate that T Tau 
Sb has been, and is still, bound to T Tau Sa, but longer-term observations will be necessary 
to demonstrate this. In any case, both T Tau Sb and Sc experienced large orbital-motion 
changes in the 1998-2001 period, when they were closest together in projection. We suggest 
that their mutual gravitational interaction at this time has resulted in these orbit changes^. 
This would amend the suggestion by Loinard, Rodriguez, & Rodriguez (2003) that T Tau 
Sa has a close companion (separation < 2 AU) whose interaction with the radio source T 
Tau Sc (which they identify as T Tau Sb) has resulted in the ejection of the latter; instead, 
T Tau Sb is a somewhat more distant companion that has probably caused the ejection of 
T Tau Sc, which is a fourth stellar or sub-stellar member in the T Tau system. 
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Fig. 1. — Narrow-band 2.26 fim (AA = 0.06 /im) image of the T Tau system, on 24 December 
2002: T Tau N (top), T Tau Sa (bottom left), and T Tau Sb (bottom right). The image 
is presented in the normal orientation (north up, east left) at a scale of 25 mas per pixel, 
and point sources are 102 mas in diameter (FWHM). Diffraction rings surround each of the 
objects; the brightest ring has a radius of about 150 mas. See Table 1 for a list of coordinates 
and magnitudes. 
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Fig. 2.— 2.145 iim (K,) image of T Tau Sa (left) and Sb (right), on 24 December 2002. 
Orientation and pixel scale are the same as in Figure 1. 
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Fig. 3. — Positions of the compact radio source T Tau Sc (solid diamonds) and the infrared 
star T Tau Sb (soHd squares), in the rest frame of T Tau Sa (sohd triangle), labelled by epoch. 
The positions for T Tau Sc are those reported for the radio source by Loinard, Rodriguez, & 
Rodriguez (2003), shifted east by 4.5 mas for reasons discussed in the text. A model orbit, 
also discussed in the text, appears as a dashed curve, with the points corresponding to the 
VLA radio observation epochs appearing as empty circles. The positions of T Tau Sb come 
from Koresko (2000), K5hler, Kasper, & Herbst (2000), Duchene, Ghez, & McCabe (2002), 
and the present work (south to north). Solid curves connecting the observations of T Tau Sc 
and of T Tau Sb are included only to guide the eye. A hypothetical epoch-2002.98 position 
for T Tau Sc, extrapolated linearly from the two most recent VLA observations, is plotted 
as an empty diamond. 



Tabic 1. Epoch 2002.98 offsets and magnitudes. 



Star Aq; (mas) Ad (mas) m(2.26 //m) 



T Tau N +20 ± 5 +691 ± 5 5.6 
T Tau Sa 8.2 

TTauSb -103 + 5 +33 + 5 8.2 



